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ABSTRACT 

O 

We investigate the high-energy photon emission around Gamma-Ray Burst 
(GRB) remnants caused by ultrahigh-energy cosmic rays (UHECRs) from the 
GRBs. We apply the results to the recent report that the supernova remnant 
W49B is a GRB remnant in our Galaxy. If this is correct, and if GRBs are sources 
of UHECRs, a natural consequence of this identification would be a detectable 
^ | TeV photon emission around the GRB remnant. The imaging of the surrounding 

emission could provide new constraints on the jet structure of the GRB. 
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Oh 1 1. Introduction 

6. 

£_! , Past Gamma-Ray Burst (GRB) occurrences in our galaxy may have left radio remnants 

resembling hyperenergetic supernovae (e.g., Perna & Loeb 2000). Recently, Chandra X-ray 
observations have led to infer that the supernova remnant W49B (G43.3-0.2) may be the 

K/ first remnant of a GRB discovered in the Milky Way 4 This remnant is located at d ~ 10 

kpc from earth, has a radius of R ~ 10 pc, and is estimated to have occurred ~ 3000 yr ago 
in our Galaxy. 

In this Letter we discuss the high energy implications of the GRB origin of some 
supernova- like remnants, and in particular of the assumption that W49B is a GRB rem- 
nant. We show that, if we also assume that GRBs are sources of ultrahigh-energy cosmic 
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rays (UHECRs), significant TeV gamma-ray emission may be detectable around W49B. A 
positive detection would verify not only that W49B is a GRB remnant, but also that the 
GRBs can efficiently accelerate UHECRs. 

GRBs are one of the most promising candidates for the origin of UHECRs (Waxman 
1995; Vietri 1995; Milgrom & Usov 1995). The physical conditions required in GRBs to 
produce MeV gamma-rays also allow protons to be accelerated up to ~ 10 20 eV. The energy 
generation rate per decade of UHECRs in the range 10 19 -10 21 eV is similar to the 7-ray 
generation rate of GRBs in the BATSE energy band 0.02-2 MeV, 

(in 
Ep—^r- (10 19 - 10 21 eV) ~ e^ RB (0.02 - 2MeV) , (1) 

and equation (1) is the common assumption made for GRBs as the sources of UHECRs 
(Waxman 2004; Vietri, De Marco, & Guetta 2003). Since photon emission is proportional 
to the square of the electron Lorentz factor, two decades in photon energy correspond to 
one decade in particle energy, but the electron (and proton) energy distribution would be 
expected to be broader than this. Also, additional photon energy is observed in many 
cases beyond the BATSE range. Thus equation (1) implying s9^ 19 10 2ii eV ~ 3eS^g 22 , MeV 
reflects approximate equipartition between nucleon and electron energies; and larger nucleon 
energies relative to electrons would not be unreasonable (e.g., Wick, Dermer, & Atoyan 
2004). Note that these energies refer to the prompt component, whereas snapshot fit and 
radio calorimetric energies refer to the afterglow. We note also that, for the purposes of this 
Letter, even UHECR of energies well below 10 20 eV can lead to GeV-TeV photon emission. 

If the GRBs are the origin of UHECRs, an unavoidable component of the UHECR 
outflow will be in the form of neutrons. This is due to conversion of protons into neutrons 
via photomeson interactions, p'-f — ► mr + on the GRB photons. The optical depth for this 
process is about a few for protons with energy > 10 16 eV (Waxman & Bahcall 1997; Vietri 
1998). Protons, unlike neutrons, feel the magnetic field in the ejecta and are subjected to 
adiabatic cooling as the ejecta expands (Rachen & Meszaros 1998), hence the neutron to 
proton ratio at similar energies could range from approximate equality to about 20% (Atoyan 
& Dermer 2001). These high energy neutrons decay into protons and electrons outside the 
ejecta, n — ► p + e~ + v e , followed by interactions with the surroundings. This results in 
a photon, neutrino and neutron emission around the remnant (Dermer 2002; Biermann, 
Medina Tanco, Engel, & Pugliese 2004). 

In this Letter we consider the photon emission from /3-decay electrons, and discuss its 
detectability in particular for the W49B case. The TeV emission is probably detectable by 
the future atmospheric Cherenkov telescope VERITAS, while HEGRA may have already 
detected this signal in their archival data, if the GRB was intense, and we strongly urge the 
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reanalysis of the HEGRA database. 



2. /3-decay electron emission around the GRB remnant 

We first consider each decade of neutron energy, to simplify the discussion. We assume 
that the GRB ejects neutrons with a Lorentz factor 7 n ~ 10 8 7 ni g and a total energy E n ~ 
10 50 i? ni 5o erg, at a time t agc ~ 3000t agc3 5 yr ago (later we integrate over the energies). These 
neutrons decay, n —>■ p + e~ + z/ e , over 

tdecay ~ Intfi ~ 3 X lO^g yr, (2) 

producing high energy electrons with a Lorentz factor 7 e ~ 7 n . Here tp ~ 900 s is the 
/9-decay time in the comoving frame. The energy of the /3-decay electrons is m p /m e ~ 10 3 
times smaller than that of the neutrons. We consider only electrons with 7 e ~ 7„ > 10 6 , 
from neutrons which decay outside the remnant, ctdccay > R, so that we can easily separate 
the /3-decay emission from the remnant emission, such as due to the pion decay (Enomoto 
et al. 2002). 

The /3-decay electrons radiate via synchrotron and the inverse Compton (IC) emission. 
The corresponding fluxes depend on the energy densities of the magnetic and photon fields. 
We adopt the typical magnetic field in our Galaxy B ~ 3/iG corresponding to a magnetic en- 
ergy density Ub ~4x 10~ 13 erg cm -3 . For the IC process, we consider the cosmic microwave 
background (CMB) as the main target photons. 5 The energy density and temperature of 
the CMB are about £/cmb ~ 4 x 10~ 13 erg cm -3 and e CMB ~ 2.7 K. The electron cooling 
time due to the combined emissions is 

t C ooi(7e) ~ t— -j- 1 x 10 4 7 - 8 1 C/ t - t 1 al ,_ 12 yr, (3) 

where C/ to tai = Ub + Uqmb = 10~ 12 £/ to tai,-i2 er S cm -3 is the total energy density, and we use 
^cmb = ^cmb min(l,m e c 2 /9ecMB7e) to take the Klein-Nishina (KN) effect into account. 6 

There are two cases, depending on whether or not most of the neutrons decay within the 
remnant age, the latter case being divided into two sub-cases depending on whether or not 
the electrons cool slower than the remnant age, leading to three different types of emission 
region shape. 



5 Our calculations show that the IR and optical backgrounds with mean values are negligible. 
6 We numerically found that the IC spectrum vF v peaks at ~ 9ecMB7e in the Thomson regime. 
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(a) Fast decay, £dccay(7e) < ^agc^ In this case the fraction of neutrons that have decayed is 
fp ~ 1. The radius of the emitting region is roughly equal to the /3-decay length ~ ct decay ~ 
0-97n,8 kpc, which is much larger than the Larmor radius of the /3-decay electrons, 

7 e m e c 2 9i 

r L ~ ^- ~ 6 x 10- 2 7e , 8J B^ pc, (4) 

qt> 

where B = 10 _6 -B_6 G. Thus, the sky distribution of the emitting region is roughly preserved, 
as shown in Figure 1 (a). The emitting region has an elongated shape on the sky, whose 
solid angle is approximately 

n ~ A9(ct dccay /d) 2 ~ (10_i7e,8)° x (10 7ei8 )° ~ 3 x lO- 3 0_ l7 2 8 sr, (5) 

if the opening half-angle of the jet is 9 = O.10_i. Note that the neutrons are beamed into a 
jet because of the relativistic beaming. 

An electron with an initial cooling time t coo i(7e) < ^agc cools down to the Lorentz factor 
determined by £ CO oi(7e) ~ ^age- Thus the current Lorentz factor of electrons is given by 

7e ' T+ 77 ~ ! x !0 8 min( 7ei8 , 3^ >3 , 5 U^ 12 ). (6) 

4£ ag eOTt/total / 

Irrespective of the initial 7 e , the order of magnitude of the current Lorentz factor ■% is the 
same in all the emitting region, so that the surface brightness on the sky depends weakly on 
the radius oc r _1 , and it is nearly homogeneous. 

The characteristic synchrotron frequency and the synchrotron peak flux at the observer 
are given by 

v m = |^ ~ 1 x 10 2 5_ 67e 2 8 eV, 
27rm e c 

vF v {v m ) ~ — — — t — - ~ 1 X 10 fpE n:50 U B ,-12l e ,8le,8 er S S Cm • ( 7 ) 

m p 7 e f/totai 4vrrf i t cool (7 e ) 
Similarly the characteristic IC frequency and the IC peak flux are 

z/cmb ~ min(7 e m e c 2 ,9e C MB7f) ~ 50min(7 ei8 ,2ecMB-37f, 8 ) TeV > 

t? I \ m e%U C MB ff3E n n - „ 2 -1 -1 -f&\ 

vF v (v CMB ) ~ — — — — - ~ 1 x 10 //3£n,50^CMB,-i27e,87 e ,8 erg s cmX8,) 

Tl^p 7e <-Aotal 47TU IcoolWeJ 

where ecMB = 10~ 3 ecMB,-3 eV is the target photon energy. The surface brightness of emission 
is about S = uF u /Q ~ 5 x 10 _9 ^l}-E ni 5 f/-i27e, 8 7er 8 er g s_1 cm -2 sr -1 . 

(b) Slow decay, tdcca y (7e) > ^a gc - In this case, almost all neutrons are still decaying. The 
radius of the emitting region is roughly ~ ct agc ~ li a ge,3.5 kpc, which is again larger than the 
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Larmor radius for our parameters. In this case, however, the shape of the emitting region 
depends on whether electrons cool within the remnant age or not. 

(bl) If t coo i(7e) > ^age> the electrons do not cool much, so the current Lorentz factor 
of the electrons is 7 e ~ 7 e . Thus the surface brightness of the emitting region is nearly 
homogeneous (see Figure 1 (bl)). The fraction of neutrons that have decayed is fp ~ 
t a ge/tdccay ~ l£age,3.57~"8 • The solid angle of the emitting region on the sky is about 

Q ~ A9(ct agc /d) 2 ~ (l^lt agc ,3. 5 )° x (10t agc , 3 .5)° ~4x 10-^-1^3.5 sr. (9) 

The synchrotron and IC emission is given by equations (7) and (8), respectively. The surface 
brightness of the emission is S — vF p /Vt ~4x 10~ 9 9ZiE nj50 U '-12^,3 .5 erg s^ 1 cm -2 sr -1 . 

(b2) If t C ooi(7e) < ^agc, almost all the electrons cool down, except for a smaller fraction of 
electrons at the jet head which have not cooled yet because the neutrons are just decaying, 
^dccay (7e) > ^agc- Although the number of these hot electrons is small, the flux from these hot 
electrons dominates that of the rest, since they have a larger energy and a shorter cooling 
time (see Figure 1 (b2)). Thus we concentrate on the hot electrons within the distance 
~ ct coo i(7 e ) ~ 37~gf/ t ~ t L al _ 12 kpc from the jet head. The current Lorentz factor of these 
electrons is -% ~ 7 e , and the fraction of neutrons that have decayed is fp ~ £ C ooiAdecay ~ 
37~g?7~ t 1 al _ 12 . The solid angle of the emitting region on the sky is 

to ~ mct agC /d)(ct cool /d) ~ (l^lt agc , 3 . 5 )° X (307^^1,-12)° 

~ 1 X 10- 2 ^lt ag e, 3 .57 e :8 1 ^ to tal,-12 *. (10) 

The synchrotron and IC emission is given by equations (7) and (8), respectively. The surface 
brightness of the emission is about S = vF u /Vt ~lx 10 _9 ^z}-E n)5 of/-i2t ag 1 e i 3.5 erg s _1 cm~ 2 
sr -1 . 

We may usually assume the emission to be isotropic in the lab-frame, since the Larmor 
radius is the minimum scale r L < min(ct coo i, ct agc ) for almost all parameters. For 7 e > 
10 10 min(2_B_ 6 U totaX _ 12 , 10 2 -B_6t agC) 3.5), however, the emission is beamed into the jet direction 
and the observed flux is reduced. 



3. Detectability 

We have numerically calculated the IC spectrum using equation (2.48) of Blumenthal & 
Gould (1970). We assume a blackbody CMB spectrum and a power-law neutron distribution 
N( , j n )d'j n cc / ~f~ p d'~f n with p = 2 for 10 7 < 7„ < 10 11 and p — 1 for 7„ < 10 7 as expected in 
GRBs (Waxman & Bahcall 1997), and use fpN{^f e < 7 e )rf7 e as the distribution of emitting 
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electrons, where fp includes the effects of the cooling and the partial neutron decay. We 
consider three cases for the normalization of JV(7 e ): (I) We set it equal to the geometrically 
corrected energy of equation (1), f 1Ql0 'y e m n c 2 N ( r y e )d'j e = 3 x 10 51 erg. (II) The same total 

energy is distributed over the entire range, j x '~f e m n c 2 N('y e )d'y e = 3 x 10 51 erg. (Ill) The 
energy of case (II) is further reduced by a factor 5, considering that only ~ 20% of accelerated 
protons are converted to neutrons (Atoyan & Dermer 2001; Dermer & Atoyan 2004). 

The emission from the /3-decay electrons is extended, hence it competes with the diffuse 
background radiation. Since W49B resides in the Galactic disk (/, b) = (43.3, —0.2), the disk 
emission is the main confusion source. We find that the IR-optical band (e.g., Bernstein, 
Freedman, & Madore 2002) and the X-ray band (Kaneda et al. 1997; Snowden et al. 
1995) are not suitable for the detection of the /5-decay emission considered here, because the 
diffuse background dominates. The MeV region (Kinzer, Purcell, & Kurfess 1999; Strong, 
Moskalenko, & Reimer 2000) could be a possible observing window, but appropriate detec- 
tors have not been developed in this band. Although INTEGRAL might be one possibility, 
the sub- MeV background may be difficult to subtract (Lebrun et al. 2004). 

On the other hand, the GeV band (Hunter et al. 1997) and the TeV band appear to be 
suitable windows for the observation of /3-decay IC emission. Strictly speaking, the diffuse 
background in the TeV band has not been detected. Only an upper limit exists (Aharonian et 
al. 2001; LeBohec et al. 2000; Amenomori et al. 2002). However the predicted background 
is below the power-law extrapolation of the upper limit (Aharonian & Atoyan 2000; Strong, 
Moskalenko, & Reimer 2004), so that it is probable that the TeV background is negligible. 

We compare now the predicted GeV- TeV flux to the sensitivities of various detectors. 
Figure 2 shows the flux of the /3-decay electrons and the sensitivities of GLAST, HEGRA, 
MAGIC 7 and VERITAS (e.g., McEnery, Moskalenko, & Ormes 2004). Note that the north- 
ern sky location of W49B makes it an unsuitable target for some ground-based detectors, 
such as CANGAROO or HESS. 

With HEGRA, MAGIC and VERITAS, whose angular resolution is ~ 0.1°, the source 
size should be taken into account. Since the search region has to be expanded, more 
background is included. The sensitivity is proportional to the inverse square of the back- 
ground, for background dominated counting statistics. Therefore the flux sensitivity of 
an atmospheric Cherenkov telescope to an extended source, intend ^ j g gj ven ^y ^extend _ 
FP° mt (fi/7r6^ ut ) 1//2 , where F° xtcnd is the sensitivity to a point source, and 6> cut ~ 0.1° is the 
angular cut in the point source analysis (Konopelko, Lucarelli, Lampeitl, & Hofmann 2002; 



7 http://hegral. mppmu.mpg.de/MAGICWeb/ 
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Lessard, Buckley, Connaughton, & Le Bohec 2001). In Figure 2 the dashed lines show the 
corrected /3-decay emission vF l/ ('D,/ir9^ ut )~ 1 ^ 2 , which should be compared with the detector 
sensitivity. 

From Figure 2 we see that VERITAS might be likely to detect the /3-decay emission. 
Even our most conservative case (III) is only a factor 2 below the VERITAS limit. MAGIC 
may also marginally detect the emission. GLAST appears unable to detect the signal, while 
the sensitivity of HEGRA is only a factor 2 above the flux of case (I). Thus, in the case of a 
stronger than average GRB (Wick, Dermer, & Atoyan 2004), the detection by HEGRA may 
not be implausible. In fact the observed total gamma-ray energy has a factor 2 dispersion 
(Bloom, Frail, & Kulkarni 2003). 

HEGRA has actually observed the region including W49B (Aharonian et al. 2001, 
2002). An upper limit on the TeV flux from W49B is given as ~ 0.14 Crab flux ~ 7 x 10~ 12 
erg cm~ 2 s _1 (Aharonian et al. 2002). However this limit constrains the flux only in a 
~ 0.1° circle centered on W49B, since the angular cut 6 cnt ~ 0.1° is applied in the analysis. 
Thus this limit is not stringent. In order to improve the sensitivity to the /5-decay emission, 
we should expand the angular cut 6 cnt so that all the /3-decay emitting region is included 
in the analysis. The probable size is ~ 0.1° x 1°, centered on W49B and with a shape 
as in Figure 1 (a), possibly along the direction in which more metals are ejected within 
W49B. Although HEGRA has ended observations, a reanalysis of the data on this region 
may constrain the /3-decay emission. 



4. Discussion 

The possibility of imaging the /3-decay emission region of a GRB remnant could open a 
novel way of constraining the structure of GRB jets. Currently jets are unresolved, but recent 
studies indicate that the jet structure is essential to an understanding the GRB phenomenon 
(Rossi, Lazzati, & Rees 2002; Zhang & Meszaros 2002; Zhang et al. 2004). Depending on 
the jet structure and the viewing angle, the same jet may be observed as different phenomena, 
such as short GRBs, long GRBs, X-ray flashes and X-ray rich GRBs (Yamazaki, Ioka, & 
Nakamura 2004, 2002; Ioka & Nakamura 2001). It is difficult to determine the jet structure 
from observations of the photon light curve and spectra during the GRB and afterglow 
phase, because the physical dimensions in these stages are much smaller than the decay 
lengths considered here, and the relativistic beaming prevents observing the whole angle 



of the jet. TeV imaging of the /3-decay emission could be a possible way 8 to determine or 
constrain the jet structure, since it provides a much longer lever arm to trace the inner jet. 

There are a number of uncertainties which could affect the conclusions. For instance, 
the photon energy density field may be higher than estimated from the diffuse backgrounds, 
due to contributions from the nearby H II region W49A. The latter emits about 10 51 Lyman 
continuum photons s _1 (Conti & Blum 2002), hence the corresponding photon energy density 
is ~ 10~ 14 (rf j 4/l kpc) -2 erg cm~ 3 , where d^ is the distance from W49A. Since the projected 
distance between W49A and W49B is ~ 40 pc, the IC emission within ~ 100 pc of W49B may 
be enhanced. Also, if the jet is directed towards us, relativistic effects such as a superluminal 
motion or differential Doppler boosting may be important. However, the barrel shape of 
W49B possibly suggests that the jet is not directed towards us. Nevertheless the forward- 
and counter-jet sizes may differ by factors of a few since it takes about the observed age for 
light to cross the system. 

To infer the jet structure precisely, the sky distribution of the electrons must be nearly 
preserved. Thus the diffusion length of electrons has to be smaller than the emitting region. 
The diffusion length r D ~ (^agc) 1 ^ 2 ~ IOO7 ' 7 tj gc35 pc within the remnant age t age is 
smaller than the emitting region for 7 e ~ 10 7 , if we use the diffusion coefficient of electrons 
k~9x 10 28 7 e7 cm 2 s _1 (Wick, Dermer, & Atoyan 2004). Electrons with 7 e < 10 7 may 
suffer diffusion, while for 7 e > 10 9 the emission is in the regime (b2) and the diffusion may 
increase the cooling region size by at most a factor 2. However we should note that the 
diffusion coefficient has large uncertainties. 

Even if W49B is not a GRB remnant, there may exist other, older GRB remnants, 
whose age may be < 10 5 yr, since the (collimation corrected) GRB rate is > 10~ 5 yr _1 
galaxy -1 . The same formalism may be used to discuss the detectability. As long as the 
electron diffusion length does not become a limiting factor, the /5-decay emission from such 
an older remnant may be detectable. 

We thank S. Razzaque, K. Mori, S. Park, J. Granot, E. Ramirez-Ruiz and the referee 
for useful comments. This work was supported in part by the Center for Gravitational Wave 
Physics under the National Science Foundation cooperative agreement PHY 01-14375 (KI, 
SK), NASA NAG5-13286, NSF AST 0098416 and the Monell Foundation. 



8 Othcr methods might include gravitational waves (Sago, Ioka, Nakamura, & Yamazaki 2004; Kobayashi 
Meszaros 2003). 
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Fig. 1. — The GeV-TeV emission region (shaded region) on the sky is shown schematically. 
(a) The /5-decay time is shorter than the remnant age £deca y (7e) < t agc , so that the radius 
of the emitting region is about the /3-decay length in equation (2). The surface brightness 
on the sky is nearly homogeneous, (bl) tdeca y (7e) > £ ag e and the initial cooling time is 
longer than the remnant age £ C ooi(7e) > ^ ag e- The radius of the emitting region is about 
~ ct age ~ liage,3.5 kpc. The surface brightness is nearly homogeneous. (b2) £deca y (7e) > £ ag e 
and t C ooi(7e) < t&ge- The radius of the emitting region is about ~ ct age ~ lt ag e,3.5 kpc. The 
jet head region, of size ~ ct coo i(7 e ), has a flux ~ £ a geAcooi times larger than the rest. The 
GRB remnant W49B, shown in the center, has a radius ~ 10 pc. 
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Fig. 2. — The flux of IC emission from the /3-decay electrons (solid lines) is compared with 
the detector sensitivity (bold long dashed lines). Three cases for the total energy are plotted: 
(I) Neutrons in the range 10 10 < 7 ra < 10 12 have the geometrically corrected energy 3 x 10 51 
erg according to equation (1). (II) The energy 3 x 10 51 erg is spread among all neutrons 
with 1 < 7 n < 10 11 . (Ill) The energy of case (II) is further reduced by a factor 5. The 
dashed lines are the flux of /3-decay emission multiplied by (n/7r^ cut ) -1 ^ 2 in order to take the 
finite source size into account, where Q is the solid angle of the emitting region on the sky 
and # cu t ~ 0.1° is the angular cut in the analysis. The sensitivities of HEGRA, MAGIC and 
VERITAS should be compared with the dashed lines. The remnant age is t agc = 3 x 10 3 yr, 
and the distance is d = 10 kpc. 



